A study of separability criteria for mixed three-qubit states 



Szilard SzalajE 
Department of Theoretical Physics, Institute of Physics, 
Budapest University of Technology and Economics, H-1111 Budapest, Budafoki ut 8, Hungary 

(Dated: February 15, 2011) 

We study the noisy GHZ-W mixture. We demonstrate some necessary but not sufficient criteria 
for different classes of separability of these states. It turns out that the partial transposition criterion 
of Peres [l|] and the criteria of Giihne and Seevinck 0] dealing with matrix elements are the strongest 
ones for different separability classes of this 2 parameter state. As a new result we determine a set 
of entangled states of positive partial transpose. 
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I. INTRODUCTION 

The notion of entanglement Q is the characteristic 
trait of quantum mechanics. It serves as a resource for 
Quantum Information Theory a relatively new field 
of science dealing with the properties, characterisation 
and applications (e^. quantum computation Ql, quan- 
tum teleportation [5|) of the nonlocal behavior of entan- 
gled quantum states. A variety of methods of Quantum 
Information Theory uses pure entangled states of a quan- 
tum system which can be easily prepared and which are 
easy to use to get nonclassical results. However, in a lab- 
oratory one can not get rid of the interaction with the en- 
vironment perfectly, thus the separable compound state 
of the system and the environment evolves into an entan- 
gled one, the prepared pure state of the system evolves 
into a mixed one. 

Generally it is a difficult question to decide whether a 
mixed state is entangled or not. A density operator 

representing the state acting on an 'H = Ti^ (Si T-L^ com- 
posite Hilbert space by definition is separable ^ when it 
can be written as a convex combination of products of lo- 
cal density operators, i.e. if there exists a decomposition 
of the form 



(1) 



where Q < Pi, ^^Pi = 1, and gf and gf are positive 
operators of trace one acting on subsystems H'^ and V.^ 
respectively. Classical correlations can give rise only to 
separable states in the sense of Eq. ([T]). 7*] 

A decomposition like the one of Eq. ([T]) is not unique, 
and it is difficult to decide whether for a given density 
operator such a decomposition exists at all. One can 
make some observations for separable pure states which 
can be extended to mixed states with the help of con- 
vex calculus. The separability criteria obtained in this 
way are necessary but not sufficient ones. (Or equiva- 
lently sufficient but not necessary criteria of entangle- 
ment.) On the other hand one can construct necessary 
and sufficient criteria using sophisticated mathematical 
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methods. [8| Unfortunately these criteria are difficult to 
use for general density matrices and only the necessary 
but not sufficient criteria can be used in practice. @ In 
this paper we calculate explicitly some of the necessary 
but not sufficient criteria of separability for a particular 
two parameter mixture of three-qubit density matrices. 
The form of these density matrices are simple enough 
to calculate explicitly the set of states for which these 
criteria hold. 

The organization of this paper is as follows. First 
of all in Section |TT] we briefly review the separability 
classes of three-qubit mixed states using the notions afc- 
separability and fc-separability. In Section IIIII we intro- 
duce our parametrized permutation-invariant family of 
three-qubit density matrices and make some observations 
about the separability class structure of permutation- 
invariant three-qubit mixed states. After having set the 
stage, in the next sections we investigate some criteria 
for separability classes. First in Section IIVI we consider 
our quantum-state as a 2 x 4 qubit-qudit system and we 
recall and use some bipartite separability criteria, namely 
the majorization and the entropy criteria related to the 
notion of mixedness of the subsystems (Sections II V Al and 
IIVBI respectively) , the partial transposition and the re- 
duction criteria which are particular cases of the positive 
map criterion (Sections IIV Cl and llVDI respectivelv). and 
the reshuffling criterion which in addition to the partial 
transposition criterion is the other one of the two in- 
dependent permutation criteria for two-partite systems 
(Section IIV E[) . As a next step in Section |V] we consider 
our quantum-state as a proper 2x2x2 three-qubit system 
and investigate some three-partite criteria for separability 
classes. In Section [VAI we recall the permutation criteria 
for permutation-invariant three-qubit case giving rise to 
another reshuffling criterion. Then we use some crite- 
ria using the expectation value of local spin-observables 
(Section IV Bp . swap operators (Section IV Cp and explicit 
expressions of matrix elements (Section lVPp . The latter 
makes it possible to determine a set of entangled states 
of positive partial transpose. In Section I VII we investi- 
gate the SLOCC classes of fully entangled states. A sum- 
mary is given in Section I VIII The explicit form of the 
corresponding matrices, some examples for permutation- 
invariant states of special separability classes and a de- 
tailed calculation of Wootters-concurrence of the corre- 
sponding two-qubit subsystems are left to Appendices [XI 
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IB] and [C] respectively. 



II. SEPARABILITY CLASSES 

A two-partite mixed state can be either separable or 
entangled, depending on the existence of a decomposition 
as given by Eq. ([1]). However, the structure of separabil- 
ity classes can be very complex even for three subsystems. 
To get the adequate generalization of Eq. ([1]) we recall the 
definitions of fc-separability and ai:-separability as given 
in Ref. [T^. 

Consider an A^-qubit system with Hilbert space H — 
(E) H'^ <E) ■ ■ ■ <S) T-L^ and denote the full set of states for 
this system as Pat. Let ak = (Si, 5*2, . . . , S'fc) denote a 
partition of the labels {1,2,..., N} into k < N disjoint 
nonempty subsets Sn- A density matrix is ak-separable, 
i.e. separable under the particular fc-partite split at, if 
and only if it can be written as a convex combination of 
product states with respect to the split ak- We denote 
the set of these states as I?^*" : 




FIG. 1. Separability classes for three qubits. The tinted sub- 
sets of the diagram can contain permutation-invariant states. 
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(2) 



where < pi, J2iPi = 1 ^^^^ ^)f" is a density opera- 
tor of the subsystem corresponding to S„ in the split 
(i.e. acting on 'H'^" = (8)aes„'W)- More generally, for a 
given k we can consider states which can be written as 
a mixture of -separable states for generally different 
a^'"* splits. These states are called k-separable states and 
denoted as 'D'^^°^: 



N 



Si-' 



(3) 



is a state of the sub- 
and in this 
,5*^'^) fc-partite splits 



in the split a^. 



can 



where < p^, J^iPi ^ and 
system corresponding to Sn'' 

case the = is['\si'\. 
be different for different i. 

Clearly I?^^"'"^^ c v'^^'^^ , so the notion of fc- 
separability gives rise to a natural hierarchic ordering 
of the states. The full set of states is Vn — V^^'^'^ 
and we call elements of v'^^'^^ \ 'D'p^^~^'^^ (i.e. the fc- 
separable but not fc -I- 1 separable states) ^^k-separable 
entangled'^ . We call the A'^-separable states fully separa- 
ble, the 2-separable states biseparable and the 1-separable 
entangled states fully entangled. 
Clearly, V"" 



N 



IS a convex set, and so is I?^ ^'^^ , because 



it is the convex hull of the union of I'^'^s for a given 
fc. Note that these definitions allow a fc-separable state 
not to be afc-separable for any particular split ak, and a 
state which is a^-separable for all ak partitions not to be 
fc -|- 1-separable. The existence of such states is somehow 
conterintuitive, but explicit examples for these states can 
be found in literature. (Using a method dealing with Un- 
extendible Product Bases Bennett et. al. have constructed 



a three-qubit state which is separable for all a2 but not 
fully separable [111] . Another three-qubit example can be 
found in Ref. [ij) 

Let us now consider the three-qubit case. (We adopt 
the notations of [To(.) For three qubits we have the par- 
titions: ai — 123, a2 = 1 — 23, a2 = 2 — 31, a2 = 3 — 12, 
a^ = 1 — 2 — 3 (here we use a simplified notation for par- 
titions usual in literature). With this, the partial separa- 
bility classes of mixed three-qubit states are as follows. 
(See also in Ref. and in Fig. [H) 

Class 3: This is the set of fully separable three-qubit 
states: P^— p ^ pi-2-3 

Classes 2.1-2.8: These are the disjoint subsets of 2- 
separable entangled states T>^~^°^ \ V^~^°'^ . Classes 2.2- 
2.8 can be obtained by the set-theoretical intersections 



of Vl'"^^ 



n 



2-31 



and V: 



f-^^. (See in Fig. [H) For ex- 
ample Class 2.8 is {V\-^^ n Pg"^^ n Vl~^^) \ Vl'^'^ 
(i.e. states that can be writen as convex combination of 
3 — 1 2-separable states and can also be written as convex 
combination of 2 — 31-separable states and can also be 
written as convex combination of 1 — 23-separable states 
but can not be written as convex combination of 1 — 2 — 3- 
separable states), Class 2.7 is {Vl^^^ n X'j^"^^) \ Pg"^^ 
(i.e. states that can be written as convex combination of 
3 — 12-separable states and can also be written as con- 
vex combination of 2 — 31-separable states but can not be 
written as convex combination of 1 — 23-separable states). 



Class 2.2 is V^^^^-^ \ (Pg U ) (i.e. states that 



can be written as convex combination of 1 — 23-separable 
states but can not be written as convex combination of 
2 — 31 or 3 — 12-separable states). On the other hand the 
union of the sets of a2-separable states is not a convex 



This 

,3-12\ 



one, it is a genuine subset of its convex hull 
defines Class 2.1 as Pg^'^^P \ (X>3~^^ U Vf-^^ U V\ 
i.e. states that are 2-separable but can not be written as 
convex combination of a2-separable states for any partic- 



3 



ular a2- However, we do not consider these states fully 
entangled since they can be mixed without the use of 
genuine three-partite entanglement. 

Class 1: This contains all the fully entangled states of 
the system: V^^"'^'^ \ V^''^^^ . 

III. A SYMMETRIC FAMILY OF MIXED 
THREE-QUBIT STATES 

Let Q be the mixture of the Greenberger-Horne- 
Zeilinger state, the W state and the maximally mixed 
three-qubit state: 

p = dil + g|GHZ)(GHZ| +w|W)(W|, (4) 
8 

where < d, w < 1 are real numbers giving rise to 
the probability distribution characterizing the mixture, 
i.e. d + g + w = 1. (In the following sections we plot the 
subsets of states for which the separability criteria hold 
on the (^-w-plane, i.e. we project the probability-simplex 
onto the d — plane. A point on this plane determines 
the third coordinate: d = 1 ^ g — w. Sometimes it is 
instructive to use the renormalized parameters d = d/8, 
g = g/2,w = w/3.) In Eq. (H)) I denotes the 8 x 8 identity 
matrix and the usual GHZ and W states are 

|GHZ> = i=(|000) + |lll)), (5a) 

|W>-i=(|001> + |010) + |100)). (5b) 

These two states are representative elements of the two 
different SLOCC-dasses of genuine-entangled three- 
qubit states. [14] The GHZ state is maximally entangled 
in the sense that its one-partite subsystems are maxi- 
mally mixed. On the other hand, its two-partite sub- 
systems are separable (having diagonal density matrix). 
The one-partite subsystems of the W-state are less mixed 
than the ones of the GHZ state, but its two-partite sub- 
systems are entangled with Wootters-concurrence 2/3. 
(See Appendix O) 

The GHZ-W mixture [d — Q line) is well studied: 
the three tangle |15|] with its convex roofs [l^, the 
Wootters-concurrences [13, [isl. the one tangle and the 
mixed-state CKW-inequality |15| were given for this mix- 
ture in the paper of Lohmayer et. al. p^ . These re- 
sults give an upper bound for values of these quantities 
on the whole simplex defined in Eq. if f{Qg.w) — 
Tahi'^^Pif{'ipi) where the minimum is taken over all de- 
composition J2iPi\'^i){i'i\ = 9g,w, and gd,g,w = dI/8 + 
(1 - d)ep/(i-d),u;/(i-d), and f(ip) = on product states 
then f{gd,g,w) < (1 - c?)/(£'g/(i-d),«,/(i-d))- 

The maximally mixed three-qubit state can be regarded 
in some sense as the "center" of the set of density ma- 
trices. On the other hand this state is sometimes called 
"white noise" because of its uniform spectrum. Mixing a 
state with white noise is the way to investigate the effect 



of environmental decoherence. [6| A noisy state is usually 
of full rank, so methods for density matrices of low rank 
(like range criterion [lOl, or finding optimal decomposi- 
tions with respect some pure-state measures) usually fail 
for such states. 

On the other hand, there are exact results for the GHZ- 
white noise mixture {w = line). In Ref. [2l| Diir and 
Cirac, using their results about a special class of GHZ- 
diagonal states have shown that g is fully separable if 
and only if < g < 1/5. Moreover, it follows from their 
observations that if the state is separable under a bipar- 
tition then it is fully separable, so Class 2.8 is empty for 
these states. In Ref. [2| Giihne and Seevinck gives neces- 
sary and sufficient condition of genuine three-partite en- 
tanglement for GHZ-diagonal states, which contain the 
noisy GHZ state: for 1/5 < g < 3/7 the state is bisepara- 
ble, yet inseparable under bipartitions, i.e. in Class 2.1, 
and for 3/7 < g < 1 the state is fully entangled. Unfor- 
tunately there are no such results for other subsets of the 
simplex given in Eq. Q. 

The noisy GHZ-W mixture given in Eq. (j?]) is clearly 
a permutation invariant one, hence the reduced density 
matrices of g are all of the same form: g^^ — g^^ — g^^ 
and g^ = g^ ~ g'^, where g^^ = Tra g, g^ = TT23 g and so 
forth. The explicit forms of these matrices are given in 
Eqs. (|X3)) and fXI)) of the Appendix. 

What can we say about the separability-classes of Sec- 
tion |TT] for permutation-invariant three-qubit states in 
general? Clearly, if a permutation-invariant state is in 
2?3^ for a particular a2, then it is in for every a2- 
So permutation-invariant states can not be in Classes 2.2- 
2.7, we have to investigate separability criteria only for 
Class 2.1, Class 2.8 and Class 3. (Fig. [H) 

(Note that the biseparability of a permutation- 
invariant state does not mean that the decomposition 
of Eq. ([3]) contains only permutation-invariant members, 
since if the latter holds then the state must be the white 
noise. To see this, write out a member of the decompo- 
sition with the help of the cr^ Pauli-matrices and Xi, yj 
real coefficients as g^ <^ g"^^ = ^{l -\- J2i ^i'^i) (8) |: (I (X" I -I- 

Ei Vj^j ® CTj") = 1(1 <^ I <^ I + I «) I + Ej %I ® 

(Jj®<yj+^ij Xiyjai®(7j®aj) which can be permutation- 
invariant if and only if a;^ = 0, j/j = 0. The reverse of 
this is that for permutation- invariant states in Classes 2.1 
and 2.8 there does not exist a decomposition as in Eq. ^ 
containing only permutation-invariant members.) 

The remaining question is wether the remaining classes 
can contain permutation-invariant states in general. 
Class 1 and Class 3 is clearly nonempty for permutation- 
invariant states, and for Classes 2.1 and 2.8 we show ex- 
plicit examples in Appendix [Bl 

If we consider the 2x2x2 three-qubit system as a 
2x4 qubit-qudit system then some well-known and easy- 
to-use two-partite separability criteria give rise to sepa- 
rability criteria for 2^3^, hence for the union of Class 
2.8 and Class 3. (This one is also a convex set since it 
is the intersection of convex ones.) First we investigate 
these criteria. 



4 



IV. TWO-PARTITE SEPARABILITY CRITERIA 

In this section we consider our system as a 2 x 4 qubit- 
qudit system (with Hilbert-spaces — and = 
7{^^) and investigate some criteria of 1 — 23-separabiHty 
which means the union of Classes 2.8 and 3. To do this 
we win need the spectra of the density matrix g given in 
Eq. Q and its marginals. (The explicit forms of these 
matrices are given in Eqs. (|Aip . (|A3|) and (|A4p of the 
Appendix.) Due to the special structure of g finding the 
eigenvalues of these matrices is not a difficult task. It 
turns out that all of the relevant eigenvalues are linear in 
the parameters g and w: 

Spect(g) = {d + 2g = {3 + 21g- 3w)/24, 
d + 3w^ (3 - 3.9 + 21w)/24, 

J= (3-35-3w)/24 6 times}, 

(6a) 

Spect(£<^3) = { 2d + 2w^{6-6g + 10w)/24, 

2d + g + w = (6 + 6.g + 2u;)/24, 

2d + g = (6 + 6.9-6u;)/24, 

2(1= (6-6.9 -6u;)/24 }, 

(6b) 

Spect(£-i) = { 4J + 5 + 2w = (12 + 4u))/24 (6c) 
4:d + g + w = (12 - 4u))/24 }. 



Here and in the following, we give expressions with and 
without ~. This is because the expressions with the quan- 
tities d, g, w can be expressive as they refer to the original 
mixing weights, on the other hand we plot in the g, w co- 
ordinates. 



A. Majorization criterion 

First of all we invoke the notion of majorization for 
probability distributions. Let p — (pi,...,p„) and q — 
(91, ... , 9„) be two probability distributions of length n. 
The definition of majorization is as follows. First we 
order p and q in non-increasing order (we denote this as 
p^). Then p is majorized by 9 by definition when the 
following inequality holds for all k: 



This is denoted by p ^ 9. (For k = n the inequality turns 
to equality since both sides of it are equal to 1. If the 
length oip and 9 diff'ers, one can add some zeroes to the 
shorter one.) 

The majorization is clearly reflexive (p -< p) and transi- 
tive (if p ^ 9 and q ^ r then p ^ r) but the antisymmetry 
{if p ^ q and q ^ p then p — q) holds only in a restricted 
manner: if p ^ 9 and q < p then p^ = q^. Hence the ma- 
jorization defines a partial order on the set of probability 
distributions up to permutations. It is clear that p q 
does not imply 9 ^ p, in other words there exist pairs of 
probability distributions which we can not compare by 
majorization. (For example let p^ := (1/2, 1/8, . . .) and 
q-^ := (1/3, 1/3, ... ), then p 7^ q and 9 / p.) This is why 
the majorization gives rise merely to a partial ordering. 

With respect to majorization the set of probability dis- 
tributions contains a greatest and a smallest element. 
One can check that all probability distribution p majorize 
the uniform distribution and p is majorized by the dis- 
tribution containing only one element: (i,^,...) -< p ~< 
(1,0,...). Employing the notion of majorization we can 
compare the amount of disorder contained in different 
probability distributions. If p ^ 9 then we can say that 
p is more disordered than 9 or equivalently 9 is more 
ordered than p, but there are pairs of distributions for 
which their rank of disorder can not be compared. 

The majorization of density matrices is defined via the 
corresponding majorization of their spectra. Let w and a 
be two density matrices, then uj ^ a hy definition when 



Spect(a;) -< Spect{a). 



(8) 



Now we can turn to the majorization criterion for 
two-partite systems. It has been found by Nielssen and 
Kempe j22f|, and it states that for a separable state the 
whole system is more disordered than any of its subsys- 
tems: 



g separable 



g ^ g^ and g -< g^ . (9) 



E 



l<k< 



(7) 



The rhs. of © can also be true for entangled states, but 
if it does not hold then the state must be entangled. 

Let us see what the majorization criterion states about 
the noisy GHZ-W mixture g given by Eq. (|4]). We can 
write out the rhs. of ([SJ explicitly using the spectra given 
by Eqs. (pa|) - ([5c| . then we have to decide whether the 
inequalities in ^ hold. For this we have to order the 
eigenvalues of the density matrices in non-increasing or- 
der. These orderings depend on the ranges of the parame- 
ters and it turns out that we have to distinguish between 
four cases. These cases are as follows: < g < ^w, 
|w <g<w,w<g< and |w < g < 1. It also turns 
out that in all these cases every inequality of ([7]) holds 
except three ones. These are as follows: 
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FIG. 2. (Color online) Majorization criterion for state Q on 
the gr-w-plane. (Red/light grey domain: g ~< and g -< g^'^ , 
blue/grey domain: g ^ g^ and g g'^'^ , dark grey domains: 
gy< g^ and g -/< g^^ .) 



where in columns (i) and {ii) are the first two inequalities 
of ([7]) (i.e. /c = 1,2) written on p ^ g^"^, and in column 
[Hi) are the first inequalities of ([7]) written on g < 
in all of the four cases. We can make the inequalities of 
PU)) expressive with the help of Fig. [2] It can be seen 
that in our case g -< g^^ implies g ^ g^, so the bigger 
subsystem (the trace map on smaller subsystem) gives 
the stronger condition. (This is not true in general. One 
can find a permutation invariant three-qubit state where 
g ^ g^ and g -< g^^ can hold independently.) 

The rhs. of ([9]) holds for states of parameters in the 
red/light grey domain of Fig. [21 so it contains Classes 2.8 
and 3. On the other hand, states of parameters in the 
blue/grey or dark grey domain are in Classes 2.1 or 1, but 
there can also be such states in the red/light grey domain. 
Moreover, the union of Classes 2.8 and 3 is a convex set 
inside the red/light grey domain. In the following we con- 
sider some other criteria in order to decrease the area of 
the red/light grey domain. In this way we can identify 
more states to be in Classes 2.1 or 1. But before this, 
we can make an interesting observation here. One can 
check that for the GHZ-white noise mixture (w — line) 
the majorization criterion g ^ g^ and g -< g^^ is nec- 
essary and sufficient for full-separability, moreover, the 
criterion g ^ g^ and g yi. g^^ is necessary and sufficient 
for Class 2.1, and the criterion g g^ and g -/< g^^ is 



necessary and sufficient for Class 1. (See Section Hill for 
summary of known exact results on the GHZ-white noise 
mixture.) Hence the condition of genuine three-partite 
entanglement is the violation of both majorization of © 
for the GHZ-white noise mixture. 



B. Entropy criterion 

The Renyi entropy of a probability distribution p is 
defined for all < a as 

H^ip) = j^lnY,P?- (11a) 
1 — a ^ — ' 

i 

For a = this is the logarithm of the number of nonzero 
PiS, known as Hartley entropy: 

Hail) hm H^[p) = In \{p, \ p, ^ 0}|. (lib) 
For a — > 1 it converges to the Shanon entropy: 

Hiip) lim H^ip) = Hip) = - V p^lnp,. (11c) 

a— )-l ^ — ^ 

i 

For — > oo it converges to the Chebyshew entropy: 

Hoo{p) lim Hcx{p) ^ -Inpmax- (Hd) 

a— vcxD 

The quantum versions of these are defined on density 
operators and can be calculated as the corresponding en- 
tropies of the spectrum. The quantum-Renyi entropy for 
all < a is 

So^ig) = lnTr(^)") - i7„(Spect(g)). (12a) 

1 — a 

The quantum- Hartley entropy is then: 

Soig) :== lim S^ig) = lnTkg = Ho{Spectig)) (12b) 

the logarithm of the rank of g. For a 1 it converges 
to the von Neumann entropy: 

Siig) : = hm Saig) = Sig) = —Trglng 

(12c) 

= iJ(Spect(£>)). 

For a — > oo it converges to the quantum- Cheby. shew en- 
tropy: 

Sooig) ■ = IhTi Saig) = — himax Spect(fp) 

(12d) 

= H^iSpectig)). 
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Now we can turn to the entropy criterion for two- 
partite density matrices (23l - [26| . This is an entropy-based 
restatement of the statement "for a separable state the 
whole system is more disordered than its subsystems" : 



g separable 

Sa{g) > Sc^ig^) and > (13) 



The rhs. of (IT5|) can also be true for entangled states, but 
if it does not hold then the state must be entangled. The 
entropy criterion follows from the majorization criterion 
since the Renyi entropies are Schur concave functions on 
the set of probability distributions. (That is if p ^ q 
then Ha{p) > Ha{q).) Therefore the entropy criterion 
can not be stronger than the majorization criterion. In 
the following we illustrate this with the state g given in 
Eq. Q for some particular choice of a. 

The rank of g, g^^ and g^ can be determined easily 
due to the simple form of the spectra in Eqs. (|6a|) - ([6c| . 
Hence the entropy criterion for Hartley entropy (jl2b[) 
can be readily examined, rkp — 8 if and only if d 7^ 0. 
The rhs. of ((T5)) holds for these states. It is true for all 
states that rkp^ = 2. On the line w — 1 — g {d — 0) 
we have to make distinction between the pure and mixed 
cases. If g = 1 (pure GHZ state) or w = 1 (pure W 
state) then rk g^^ — 2 and rk p = 1 hence for this case 
So{g) ^ So[g^) and So{g) ^ Soig^^). For the genuine 
mixtures of GHZ and W states rk g^^ = 3 and rk p = 2 
hence So{g) > So{g^) but So{g) ^ So{g'^^). So we can 
establish that the entropy criterion in the limit a 
(quantum-Hartley entropy) is too weak, it identifies only 
the GHZ- W mixture to be entangled. 

Consider now the entropy criterion in the a — ?■ 00 limit. 
This can easily be done because the inequalities of the 
rhs. of (|13p are the same as the ones in the (i)th and 
(m)th column of (jlOl) . which are written on the maximal 
eigenvalues. Hence in this case we have fewer restrictions, 
and one can see in Fig. [3] that the rhs. of ([T3| holds 
for more states than the rhs. of © in the case of the 
majorization criterion. Hence the entropy criterion in the 
a — 00 limit (quantum-Chebyshew entropy) identifies a 
little bit fewer state to be entangled than the majorization 
criterion. 

Increasing a from to c» one can see in Fig. 2] how 
the borderlines of the domains of the entropy criterion 
shrink to the ones in Fig. |3l It is not true in general 
that if Haip) < Haiq) and a > /3 then Hpip) < Hp[q). 
For these particular spectra it seems that the domains of 
smaller as would contain the domains of larger as, but 
for the large values of a one can see that this is not true. 
However, no line can cross the border of the domain of 
majorization criterion, since the entropy criterion can not 
be stronger than the majorization criterion. 
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FIG. 3. (Color online) Entropy criterion in the limit a — s- 
00 for state (Q on the g-iu- plane. (Red/light grey domain: 
Soo{q) > Sooig^) and Sac{g) > Soo{g'^''^), blue/grey domain: 
Sco{g) > Sao{g^) and S'oo(e) ^ S^oig^^), dark grey domains: 
S^{g) t Sooig^) and S^ig) t Soo(f?'').) 
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FIG. 4. (Color online) Entropy criterion for a = 
i, |, |, 1, 2, 3, 4, 5, 10, 20 for state (01) on the p-w-plane. 
(Red/grey curves: the border of the domain Sa{g) > ^^(^i^) 
and Sa{g) > Sa{g^^), blue/black curves: the border of the 
domain Sa{g) > Sa{g^) and Sa{g) ^ Sa{g^^). Red/light grey 
domain: copied from Fig. [3] of the a — > 00 case.) 

C. Partial transposition criterion 

If a two-partite state is separable then the partial 
transposition on subsystem A acts on the gfs of the de- 
composition given in Eq. ([T]). The transposition does 
not change the eigenvalues of a self-adjoint matrix, so 
{gf)'^s are also density matrices (i.e. self-adjoint matrices 
of trace one). Hence the partial transpose of a separable 
density matrix is also a density matrix. (Its eigenvalues 
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are not the same in general as the ones of the original ma- 
trix, but they are also nonnegative and sum up to one.) 
The reverse is not true unless the system is of qubit-qubit 
or qubit-qutrit 8] , so our 2x4 system is the smallest one 
for which this implication is only one-way. In general we 
get the partial transposition criterion of Peres [l| : 

g separable => g^^ > 0. (14) 

It is clear that no matter which subsystem is transposed. 

(The partial transposition criterion is the consequence 
of the positive maps criterion: Q 

Q separable ■4==^> 

($ (g) 1)^1 > for ah positive maps (15) 

This is a necessary and sufficient criterion, but we can not 
check it for all $. But we can consider a particular class 
of positive maps to obtain necessary but not sufficient 
criteria. For example for $(a;) = ur^ we get back the 
partial transposition criterion.) 

Let us apply the partial transposition criterion to the 
state Q of Eq. The spectrum of g^^ can easily be 
calculated due to its block-structure. (See in Eq. (jA2[) of 
the Appendix.) 
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FIG. 5. (Color online) Partial transposition and reduction 
criteria for state (j4]) on the g-ro-plane. (Eqs. (|17ap and (|17b|) 
of partial transposition criterion hold inside the intersection 
of the red/grey ellipses. Blue/black curves: the borders of do- 
mains inside the additional Eqs. (|20ap and (|20bp of reduction 
criterion hold. Red/light grey domain: copied from Fig. [2] of 
majorization criterion.) 



Spect(£<^i) = 



= (3 - 35 + w ± 4 VV+^)/24, 



d + 5/2±v/,g2 + 8^2/2 

= (3 + 3.9 - 3w ± 2v/V + 32w2)/24, 
d + g = (3 + 9.9- 3w)/24, 
d + 2w=(;i- 3.9 + 13w)/24, 

(1= (3 - 35 - 3w)/24 2 times}. 

(16) 

Only the minus-version of the first two eigenvalues can 
be less than zero hence we get two inequalities for the 
positivity of g^'^ : 



0< d^ +dw~g^ 

2 2 (17a) 

< -135g^ - 15w^ - 6gw - 18g + 6w + 9, 



r,2 



0<d^+dg- 2w 
< -27.9^ - llOw^ - 18gw 



(17b) 



Each inequality of these holds inside an ellipse. These 
ellipses intersect nontrivially and in the intersection the 
rhs. of (ITi|) holds. (Red/light grey curves in Fig. [S]) 
5 = 1/5 and w = (24\/2 - 9)/119 = 0.209589 . . . arc the 
bounds for the union of Class 2.8 and 3 for the GHZ-white 



noise {w = 0) and the W-white noise {g = 0) mixtures 
respectively. 

The partial transposition criterion states then that if 
a state is in Classes 2.8 or 3 then its parameters are in- 
side the intersection of the ellipses, but there can also 
be states of Classes 2.1 or 1 in there. On the other 
hand the states must be in Classes 2.1 or 1 for param- 
eters outside. The inequalities of |_?ya[ ) and |j7&p are 
strong in detection of GHZ and W state respectively. In 
Fig. [S] we have also plotted the corresponding domain of 
the majorization criterion. (One can check that the only 
intersection-points of the borderlines of the correspond- 
ing domains of the two criteria are (g = 2/13, w = 3/13) 
and (5 = 1/5, w ~ 0). This criterion is also a necessary 
and sufficient one for the full separability of the w = 
GHZ-white noise mixture.) It can be seen that the par- 
tial transposition criterion gives stronger condition than 
the majorization criterion, it identifies more state to be 
in Classes 2.1 or 1. Hence the majorization criterion can 
not identify entangled states of positive partial transpose 
(PPTES) on the simplex defined in Eq. (g]). 

The PPTESs are exotic entangled states. They are 
bound entangled (undistillable) states, pT j that is entan- 
gled states from which no entanglement can be distilled 
at all. (The entanglement distillation [2^ [2^ is a family 
of methods which allow one to extract locally maximally 
entangled pure states out of a given state or its copies.) 
It is usually hard to check that a state of positive partial 
transpose is not separable, there are few explicit exam- 
ples of PPTES in the literature (see a list of references 
in Section 1.2.4 of Ref. 29]). All the states in Class 2.8 
are PPTESs. 
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D. Reduction criterion 

The next one of the examined criteria is the reduction 
criterion. [30, It states that 



separable 



g > and 



~ g>0. (18) 



This is the consequence of the positive maps criterion 
given in (jl5p for the particular positive map $(w) = 
Tr(a;)I — uj. The importance of this criterion is that 
its violation is sufficient criterion of distillability. (30| 
It is known (30l | that the reduction criterion can not be 
stronger than the partial transposition criterion and they 
are equivalent for qubit-qudit systems. Since our state 
g defined in Eq. (|4]) is the permutation invariant one of 
three qubits considered as a 2 x 4 qubit-qudit system, the 
equivalence of these two criteria means that some kind of 
pure state entanglement between 1 and 23 can be distilled 
out from every state of non positive partial transpose. In 
other words in the simplex defined by Eq. ^ there are 
no bound entangled 2x4 states of non positive partial 
transpose. 

We can illustrate the equivalence of the partial trans- 
position and reduction criteria. To do this we have to 
examine the positivity of the matrices (g) g^^ — g and 
pi (g)l23 - g. (See in Eqs. (IX5|) and (IM)) of the Appendix.) 
Since Tr(a;)I — cj = (cr2W0'2)"^ for 2x2 matrices (with the 
Pauli matrix = [ ^ V ] ) turns out that 



Spect(li (g) £-23 _ : 
Spect(ei ® - g) - 
{3d-H3w/2±V4^ 



Spect(£i^i), 



(19a) 



= (9 - 9.9 + 3w ± 4V9g2 -f w2)/24, 
3d + .9 ± %/2w = (9 + 3.9 - 9w ± %V2w) /24, 
3d + g + 2w (9 + 3.9 + 7w)/24 2 times, 
M + g + w^[9 + 2,g-w)/2'i 2 times }. 



(19b) 

For P (g) g^^ — g > we have the same conditions as 
in Eqs. (|17a|) - (jl7b|) of the partial transposition criterion. 
The additional inequalities arise from the minus- version 
of the first two eigenvalues of g^ (g) I^'^ — g: 



Oil 

3 



< 9d'' - 
< -63.9^ 



7w^ ~ 54.9W - 1625 -I- 54it; 



< 3d + g- V2w 
0< 35- (9 + 8V2)w- 



81, 
(20a) 

(20b) 



The first one of them is true outside a hyperbola, the 
second one is true under a line. (Blue/black curves in 
Fig. El) 



It can be seen that the last two inequalities (I20a[) - (l20bp 
do not restrict the ones in Eqs. (|17al) - (jl7bp . as it has to 
be, and because of (|19ap the reduction criterion and the 
partial transposition criterion hold for the same states 
of the GHZ- W-white noise mixture. Here we get the 
stronger condition for the map $(0;) = Tr(a;)I — cj acting 
on the smaller subsystem. We can also observe that the 
inequalities of (|20ap and (|20b|) are good in detection of 
GHZ and W state respectively, but not so good as the 
ones of partial transposition criterion. However, one can 
check that on the w = GHZ-white noise mixture the 
reduction criterion (g) g^^ — g> and g^ ^t^^ — g > is 
necessary and sufRcient for full-separability, the criterion 

(g) g'^^ — g ^ and (g P'^ — g > is necessary and 
sufRcient for Class 2.1, and the criterion (g g^^ — g ^ 
and g^ (E) — g ^ is necessary and sufficient for Class 
1 in the same fashion as in the majorization criterion of 
Section IVAl 



E. Reshuffling criterion 

The reshuffling criterion is independent of the partial 
transposition criterion, so it can detect entangled states 
of positive partial transpose. It states that: 



g separable 



|i?(^?)l|Tr < 1, 



(21) 



where the trace-norm is PIItt = TrvMU, and the 
reshuffling map R is defined on matrix elements as 

[R{g^]ii' Jj' Qij^i'j'- 

The four nonzero singular values of the 4x16 reshuffled 
density matrix, (see in Eq. \A7\ of the Appendix) i.e. the 
square root of the nonnegative eigenvalues of R{g)R{g)'^ 
are: 



Spect ^R{g)R{gy 
= l^pi ± 2^/2, + 2w^, + 2w^Y (22) 



where 



Pi =l6d'^ 

P2 =64d'* 



4.9^- 
9w^ 



lOw^ + 8dg + 12dw, 
- Q4d^g + 96(Pw + 12dw^ 



Wd^g^ 



80(P gw -I- Iddg'^w -f 2Adgw'^ . 



The sum of them is less or equal than 1 inside a curve of 
high degree which can be seen in Fig. |6] (red/grey curve). 
States of Classes 2.8 and 3 must be inside this curve, 
states outside this curve must belong to Classes 2.1 or 1, 
but one can see that this criterion does not restrict the 
partial transposition criterion, it does not detect PPTESs 
in the GHZ- W-white noise mixture of Eq. (U). 



9 




0.2 



0.4 



0.6 



FIG. 6. (Color online) Reshuffling criteria for state Q on the 
g-w-plane. (Red/grey curve: reshuffling criterion for 2x4 
system as in Section FlV El green/light grey curve: reshuffling 
criterion for 2x2x2 system as in Section IV Al We have 
also copied the borderlines of the domains in which (|17a|l and 
()17b|) of the partial transposition criterion hold from Fig. \5\ 
The inequalities hold on the side of the curves containing the 
origin.) 



V. THREE-PARTITE SEPARABILITY 
CRITERIA 

In this section we consider our system given in Eq. ^ 
as a proper 2x2x2 three-qubit one and investigate some 
general 3-qubit /c-separability criteria. 



A. Permutation criterion 

First consider the permutation criterion in general 
[33 |. Note that the reshuffling and the partial transpose 
of a density matrix are nothing else than the permu- 
tation of the local matrix indices. Moreover, since 
the trace norm is the sum of the absolute values of 
the eigenvalues for hermitian matrices and the trace is 
invariant under partial transposition it turns out that 
£1^1 > if and only if |j£i^i||Tr = 1- So the partial 
transposition criterion and the reshuffling criterion 
mi can be formulated in the same fashion. Moreover, 



this can be done for N subsystems in a general way. [32 1 
Let TT e S'2Ar a permutation of the 2N matrix indices and 
let the map realizing this index permutation: ii g — 

S Qili2...iN,iN + liN + 2---i2N 1*1*2 • • • iN){iN+liN+2 ■ ■ ■ i2N\ 

then [^ir(f?)] V(i)i^(2) ...i,r(JV) ,i„(JV-|-l)i„(JV + 2)---i,r(2W) ~ 

ftii2...ijv,ijv+iiN+2...i2jv Now the permutation crite- 
rion states that 



Q fully separable 



|A.(e)||Tr < 1, V^e52iV. (23) 



The permutation criterion gives {^N) ! criteria but not 
all of them are inequivalent. It is well known (32i] that for 
two subsystems, every criteria given by the permutation 
criterion turn out to he equivalent either the partial trans- 
position criterion or the reshuffling criterion. In Ref. ^33j 
Clarisse has shown that there are only six inequivalent 
criteria in the case of three subsystems: three one-partite- 
transpositions and three two-partite-reshufflings. For our 
permutation-invariant three-qubit system all the one- 
partite-transpositions give the same condition which we 
have already investigated in Section IIV CI On the other 
hand, all the two-partite-reshufflings give another condi- 
tion which is a new one. 

So let R' = the map implementing the reshuffling 
of the 2 and 3 subsystems: [R' {Q)]ijj\i'kk' = Qijk,i'j'k'- 
(See in Eq. (jA8|) of the Appendix.) With this we have 
to calculate the eigenvalues of the matrix R' {g)R' {qY 
for the two-parameter state g given in Eq. (jj]). This 
8x8 matrix can be transformed (by simultaneous row- 
column permutation) into blockdiagonal form consisting 
three blocks of the shape 3x3,3x3 and 2x2. However, 
the forms of the w-depending eigenvalues of the 3 x 
3 blocks are still too complicated, so we only plot the 
border of the domain in which the criterion holds. 
(See green/light grey curve in Fig. [SI) 

The condition ||i?'(£))||Tr < 1 holds inside the 
green/light grey curve in Fig. [5| This figure suggests 
that this reshuffling does not give stricter condition for 
full separability than the partial transposition criterion, 
hence it can not identify PPTESs. However, we can 
not be sure in this due to the difficult computation of 
||i?'(g)||Tr- Fully separable states must be enclosed by 
the curves belonging to (|17ap - (|17bl) of partial transpo- 
sition criterion, states outside this domain must belong 
to Classes 2.8, 2.1 or 1. However, in Section [IV CI the 
partial transposition criterion has yielded condition for 
Classes 2.8 and 3, so we can conclude that states outside 
this domain must belong to Classes 2.1 or 1, the Class 
2.8 is completely restricted into this domain. 



B. Criteria on spin-observables 

In Ref. [l3] Seevinck and Uffink introduced a system- 
atic way to obtain necessary criteria of separability for 
all the separability-classes of an A^-qubit system. Their 
criteria generalize some previously known criteria, (see 
Refs. in Ref. [10]) such as Laskowski-Zukowski crite- 
rion (necessary for fc-separability), Mermin-type separa- 
bility inequalities (necessary for fc-separability) , Fidelity- 
criterion (necessary for 2-separability) and Diir-Cirac de- 
polarization criterion (necessary for ckfe-separability). We 
consider the three-qubit case and get criteria for Class 
2.1, Class 2.8 and Class 3 given in Section [ill 

The method of Seevinck and Uffink deals with 
three orthogonal spin-observables on each subsystem: 
(X(i),r(i),Z(i)). Here the superscript (1) denotes that 
these are single-qubit operators. Let I^^^ denote the 
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2x2 identity matrix. From the (X^^), F^^), Z^^', /(i)) 
one-qubit observables acting on the 2 and 3 subsys- 
tem one can form two sets of two-qubit observables: 
(xF',yi^\zi^\/l^'). Here the superscript (2) denotes 
that these are two-qubit operators and a; = 0, 1 refers to 
the two sets: 



X 

Z- 
I 



(2) 



(2) 



(2) 



(2) 



-(2) = 1 (^/(i) ^ Jd) + ^(1) ^ ^(D) 



(24) 



(Note that Ix s are noi identity operators.) From 
this two-qubit observables and the one-qubit ones 
acting on the 1 subsystem one can form four sets 
of three-qubit observables acting on the full system: 
{x'^\Yi'^\ z'i'\l'^^). Here the superscript (3) denotes 



that these are three-qubit operators and x 
refers to the four sets: 



y 



X. 



(3) 

y+1 



y(3) 

y 



Y, 



(3) 



7(3) _ - 

y 

A3) _ " 



y(i)( 
y(i)( 



X 



X 



(2) 
a/2 

(2) 
a/2 



(2) 



y/2 



(2) 



y/2 



>Y 



(2) 



a/2 



(2) 

a/2 



r(2) 

a/2 

r(2) 

a/2 

7-(2) 
y/2 " 



(2) 



a/2 



J//2 



7(2) 

"a/2 

7(2) 

"a/2 

7(2) 
"a/2 

7(2) 
"a/2 



(25) 



for y = 0, 2. (Note that Ix s are not identity operators.) 

Now for particular a2 investigating some relations 
among the expectation-values of these operators with re- 
spect to the state g one can get some nontrivial inequal- 
ities valid for all g G T)^^ . From these one can form 
inequalities valid for a given separability class of Section 
im Here we recall IcI these criteria for the classes we 



need to deal with: 

^T^2— sop 

g ePg 



+ (ri'V < E - (26) 

y^x 



for X = 0,1,2,3, 



02 

3 



nax{(X(3))2 + (yj3)^2| 

<mm{(43))2-(Z(3)>2}<l/4 



(27) 



and 



_^3— scp 

g gPg 



c{(x(3))2+(rj3)^2| 

<mm{(/(3))2-(Z(3))2}<i/i6. 



(28) 



One has to do some optimisation of the local spin observ- 
ables {X'^^\Y^^\ Z'-^^) to get violation of the respective 
inequality for a given state. 

In the following we will consider some spe- 
cial measurement-settings when the observables 
are the same for each subsystem. 

Writing out explicitly (Xx^^ , Yx^^ , zi^^ , li^^ ) , one can see 
that for a permutation-invariant state the squares of 
the expectation values are the same for x = 1,2,3, i.e.: 



0,1,2,3 {X[ 



(3)\2 



{X: 



(3)\2 



{X^^^y, and the same for Yx'^'s 



(3)c 



Z'i^-'s and li^^s. Hence we have to consider merely the 
x = 0, 1 indices. 

First consider the special choice when 
(X(i),r(i),z(i)) = iax,ay,a^) for each subsystem. 
(Here CiS denote the usual Pauli matrices.) We refer 
to this as Setting I. The inequalities can be 

written as relatively simple expressions in the matrix 
elements uM: 



g eP'-'^'P 



IpotI < y/geegu + V£'556'22 + y/gssgu, 



(29) 



leeil < y/eaogn + Vfi'556'22 + y/gssga, 
\g52\ < y/geegu + y/googri + y/gsagu, 
\e3i\ < VeeePii + Vf'556'22 + y/gmgrr, 



gef]vr 



a2 

max 



{l^'07|M^?6l|M^>52|M^^34|'} 



< min {goQg77, geegu, g55g22, £'33 6'44} < 1/16 



(30) 



and 



_ ^3— scp 

gev^ 



maxjleoyP, IffeiP, |£'52p, |£'34p} 

<mm{googr7, geegii, g55g22, g33g44} < 1/64. 



(31) 
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(Here the matrix indices run from to 7 so as to be 
equal to the binary indices we use later.) Let us con- 
sider another two special measurement settings: Setting 
II: = {ay,a^,a^) for each subsystem, 

Setting III: F^^), Z^^)) = {(T^,a^,ay) for each sub- 

system. The inequahties of (|26)) - (|28| written for these 
two settings are much more complicated expressions in 
symbolic matrix elements than the ones in (P9|) -([5T |) . But 
for the state g given in Eq. it is not too difficouh to 
write out these inequalities explicitly. It turns out that 
for each of these three settings the x — 1 inequality of 
(|26l) . the second inequality of ([27|) and the second in- 
equality of (|28)) hold for all the parameter values of the 
simplex. Because of this, the criteria hold for Class 3 
are not stricter than the ones for the union of Class 2.8 
and Class 3. The remaining inequalities for the three 
measurement settings are as follows: 



^ T^2— sen 



II. 
III. 

and 



g < 3\Jd{d + w) 

< -7g^ - 6gw - 15w^ - 18g + 6w + 9, 
3w < \/{8d + w){8d + 4:g + w) 



^J^~g^ + 81ii)2 < 3(8J + 2g + w) 



< -g^ - hw^ - 12u;-|-9, 



(32a) 

(32b) 
(32c) 



,a2 



11. 



111. 



< d{d + w) 
< -45.g2 - 2gw - 5^^ - 65 + 2w + 3, 
81w^ < {8d + w){8d + Ag + w) 

< -9.9^ - 77^2 - 12w + 9, 
4g2 + 81w^ < {8d + 2g + wf 

< -9.g^ - 77w2 - \2w -f 9, 



(33a) 
(33b) 
(33c) 



Clearly, the inequality of (|32cl) is weaker than the one of 
(j32bp . the inequality of p3cp is the same as the one of 
(j33bp . Moreover, the inequality of p3ap is the same as 
the one of (|17ap of partial transposition criterion, but the 
inequality of p3bp is strictly weaker than the other one 
of partial transposition criterion. So these settings does 
not give stricter conditions for Classes 2.8 and 3 than 
partial transposition criterion, however, we get criteria 
for biseparability for the first time. In Fig. [7] we show 
the borderlines of the domains of the criteria belonging 
to Settings I. and II. These inequalities restrict Classes 
2.1, 2.8 and 3 to be inside the domain enclosed by the 
blue/black curves and Classes 2.8 and 3 to be inside the 
domain enclosed by the red/grey curves. We can con- 
clude that Settings I. and II. are strong in detection of 
GHZ and W state respectively. One can check that for 
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FIG. 7. (Color online) Criteria on spin-observables for state 
(|4]) on the g-w-plane. (Red/grey curves: the border of do- 
mains inside Eqs. (|32a|) and (|32b|) hold, blue/black curves: 
the border of domains inside Eqs. (|33ap and (j33b[) hold. The 
inequalities hold on the side of the curves containing the ori- 
gin.) 



the w = GHZ-white noise mixture the inequalities of 
P2ap and p3ap of Setting I. hold if and only if the state 
is fully separable, (I32ap is violated but (I33ap holds if and 
only if the state is in Class 2.1 and both of them are 
violated if and only if the state is fully entangled. For 
the g = W- white noise mixture if 3/11 < w then g is 
in Class 2.1 or Class 1, and if 3/5 < w then g is fully 
entangled. 

However, there are infinitely many criteria depend- 
ing on the measurement settings and we do not have a 
method to find a set of settings leading to the strichtest 
criterion. We have tried some other randomly chosen set- 
tings which can be used to reduce the area where the cri- 
teria hold. We could not find settings that give stronger 
criteria on the w = oi g — axes of the simplex than 
Settings I. and II. respectively. We have found settings 
that excludes states from the corresponding classes, but 
these states are far from these axes, and we have not 
found settings which give stronger condition for Classes 
2.8 and 3 than the partial transposition criterion. We 
have found settings by which the condition for bisepara- 
bility can be strengthened, but these conditions are just 
a little bit stronger far from the axes than the ones in 
Section |Vni 



C. Criteria on matrix elements 

In a recent paper |34j, Gabriel et. al. have given cri- 
terion for fc-separability, based on their previously de- 
rived framework for the detection of biseparability [s^ . 
It turns out that for the noisy GHZ-W mixture given in 
Eq. ^ these criteria give the same results as the ones of 
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Seevinck and Uffink, given in the previous section, but 
these criteria have the advantage that they can be used in 
the same form not only for qubits, but for subsystems of 
arbitrary, even different dimensions. As our knowledge, 
these are the only such criteria of fc-separability. 

Consider some permutation operators acting on 
H (E) H, i.e. on the two copies of the A^-partite 
Hilbert space V. = rO- ® H"^ ® ■ ■ ■ ® . Let 
PaS the operators which swap the ath subsystems 
of the two copies: Pa\iii2 ■ ■ -iN) ® \j1j2 ■ ■ ■ Jn) = 

\ili2 ■ ■ -ia-ljaia+l ■ ■ -iN) "Si |jlj2 ■ • • Ja-l^aja+l • ■ • jjv) 

where {|ia)} is basis in T-L"". Now for a composite 
subsystem = (^a^sH" let Ps = ]\aesP<^- The key 
fact is that if the state of a set of subsystems can be 
separated from the rest of the state then the corre- 
sponding Ps leaves the two copies of the state invariant: 
PgQ®'^Ps = Q®"^ ■ With this and convexity arguments, 
one can get the following criteria for fc-separability [33 |: 

^ ^fc — sen 

Ck \ A 

ri=l " / 

(34) 

where |$) G 'H 8 'H is a fully separable vector, and the 
total swap operator is Ptot = Il^i ^a- Here i runs over 
all posible fc-partite splits a\^'' = {s["\si"\ S^!^). 

The inequality in p4p is written on the matrix ele- 
ments of g determined by the separable detection- vector 
1$). For a given state, optimisation on |$) is needed to 
achieve the violation of (|34p . 

To apply these criteria to the noisy GHZ-W mixture 
given in Eq. (jH) we have to choose a suitable detection- 
vector |<i>). It turns out that |$ghz) = |000111) and 
|$w) = H'^^l^Quz) are good choices for states in the 
vicinity of GHZ and W states respectively, as observed in 
Rcf. [3^. (Here H = 1/^2 [ } _\ ] is matrix of the usual 
Hadamard/discrete Fourier transformation for qubits.) 
With these two vectors we get the same criteria for 2- 
separability as the ones in (I32a[) and (I32b[) respectively, 
and for 3-separability as the ones in p3al) and (j33b[) re- 
spectively. However, (|33a|) and ()33bl) obtained by the cri- 
teria on spin observables are criteria not only for Class 
3, but for the union of Classes 2.8 and 3, so in this sense 
the criteria on spin observables are a bit stronger. 

We can not be sure that the detection-vectors above 
give the strongest conditions at least for the noisy GHZ 
and noisy W states. However, it is an interesting ob- 
servation that the Hadamard transformation relates not 
only the two "optimal" detection- vectors I^ghz) and 
|$w) but also the two "optimal" measurement-settings 
Setting I. ((Ta;, (Ty, (Tz) and {a^, —<Jy,<Jx), (by the trans- 
formation ai f-7> Ha.iH^) which is equivalent to Set- 
ting II. (This equivalence holds only for permutation- 
invariant three-qubit states, when the three sets of ob- 
servables (A:(i),r(i),z(i)) are the same for each sub- 
system. In this case one can check that the quantities 



{XfY + (yi^V for a: = 0,1,2,3 are invariant under 
the transformation {X^^\Y''^\ Z^^^) o {Y^^\X^^\ Z'^^^) 
and (A:(i),y(i),Z(i)) o (A:(i),-r(i),Z(i)). These can 
be seen by writing out the definitions given in Eqs. (j25p .) 

We have tried some other randomly chosen detection- 
vectors which can be used to reduce the area where the 
criteria hold, and we get the same observations as at the 
end of the previous section: One can strenghtend the 
conditions only far from the = or 5 = axes of the 
simplex, we have not found detection-vectors which give 
stronger condition for full-separability than the partial 
transposition criterion, and we have found settings by 
which the condition for biseparability can be strength- 
ened, but these conditions are just a little bit stronger 
far from the axes than the ones in Section IV Dl 

D. Criteria on matrix elements — a different 
approach 

In Ref. [2I] Giihne and Seevinck have given some further 
biseparability and full-separability criteria on the matrix 
elements: 

^ T^2 — sen . 

(35a) 

|£'12| + Igiil + \g2i\ < \/ 900 033 + y/QmQ55 + y/m066 

+ {gii + Q22 + £'44) /2. 

(35b) 

The criterion p5ap is necessary and sufficient for GHZ- 
diagonal states and can also be obtained as a special case 
of the criteria of Section FV Bl (Eq. ([29])). However, this 
criterion — and the others in this Section — arises from a 
quite different approach as the one in ([29]), since these cri- 
teria have been derived from direct investigation of the 
matrix elements of pure separable states with the use of 
convexity argument. The criterion in (|35b|) is indepen- 
dent of the first one and it is quite strong in detection of 
W state mixed with white noise. 

Of course these and the following inequalities can be 
written on local unitary transformed density matrices, 
and optimisation with regard to local unitaries might be 
necessary, but this can lead to very complicated expres- 
sions in the original matrix elements. An advantage of 
the method of the previous section is that it handles the 
matrix indices through the detection vector | $) . 

The full-separability criteria of Q : 

^ scp 

IQOtI < (£'ll£'22e33£'44e55£'66)^^'' , (36a) 
|£'12| + |eil4| + |£'24| < V £'00^33 + V6'00£'55 + 

(36b) 

P6ap is necessary and sufficient for GHZ state mixed 
with white noise, and (|36bp is violated in the vicinity of 
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the W state. Moreover, one can obtain other conditions 
from p6al) by making substitutions as foUows. Consider 
a fuUy separable pure state g^-scp _ |v]/) (\l>|^ where l^*) = 
(ao|0) + ai|l)) ® (6o|0) + 6i|l)) (co|0) + ci|l)). Then 
the diagonal elements are {Q^~^''^)ijk,ijk = |aiPl^jP|cfcP, 
where we use the ijk = 000, 001, . . . , 111 binary indexing. 
Then 

/ 3— scp\ / 3— sop~i 

\t: )ijk,ijk\f^ }i' j' k' ,i' j' k' 

('/,3-sopN (' 3— sep\ 

\t! )ijk' ,ijk' \b! k^i j k 

( „3 — scp\ ('^3— scp\ 

k.ij' k\Lf )i'jk',i'jk' 

— \Lf )i'jk,i'jk\y k' k' ■ 

Moreover, the rhs. of the inequality of p6al) can be writ- 
ten as (f?ii£'22£'33£'44£'55£'66) ^ ^^-^ with these substi- 
tutions we can obtain a third power of four matrix ele- 
ments under the l/12th power. So we can get expres- 
sions of four matrix elements (e.g. {QiiQ22QAAQn)^^^) on 
the rhs. of (I36ap . With the substitutions above one can 
get 28 different inequalities for (j36ap with an expresi- 
son of sixth order under the sixth root on the rhs. and 
12 different ones with an expresison of fourth order un- 
der the fourth root. For permutation-invariant states we 
have Qii = Q22 = £'44 and gee = Pas = ^33 and for our 
case QoQ — Qn as well, hence the number of different 
inequalities reduces to 8 and 5 respectively. The rhs's 
of inequality p6a|) which are different for permutation- 
invariant matrices are as follows: 

UooeTTj , (£'ll£'22£'33£'44£'55£'66) 
/ 2 \ 1/6 / 2 \ 1/6 

\QmQ\\Q22Q^'^e^^) , [Q00933Q55Q66Q77) , 

( 22 \ 1/6/2 2\l/6 

(£'00£'ll£'66£'77j , (£'oo£'ll£'66£'77j j 

12 ^l/6 / 2 \l/6 

\QllQ22QiiQmQnj , (,£'00£'ll£'33£'55£'66j 

and 

( ^22 6133 ^44 ess) , (£'00£'lie'66£'77j , 
(£'ll6'22£'44 6'77j , (£'00£'336'55£'66j • 

It turns out that the strongest conditions can be given 
with the last of these and with the original one in p6ap . 
(We could also make some substitutions in the rhs. of 
(j36b[) but these would not give stronger conditions than 
the original one.) 

Writing out the criteria of biseparability and full sep- 
arability we get: 

g e P^-scp ^ ~ ^ 3^Jd{d + w), (37a) 
w < \/{d + g)d+{d + w)/2 (37b) 




0^2 04 0^6 0.8 i 

g 



FIG. 8. (Color online) Criteria on matrix elements for state 
@ on the g-m-plane. (Green/light grey curves: the borders of 
domains inside Eqs. (|38a|) - (|38c|) hold, blue/black curves: the 
borders of domains inside Eqs. (|37al) - (l37b|) hold. Red/grey 
curve: the border of domain inside Eq. (|17b|l of partial trans- 
position criterion hold, copied from Fig. [5] The point g = 1/5, 
w = 1/5 is also shown. The inequalities hold on the side of 
the curves containing the origin.) 



and 

^g^3-scp ^ g<{{d + w)d)^/^, (38a) 

~g<{{d+~g)d^)'/\ (38b) 

w < \J{d + g)d. (38c) 

(See in Eq. (jAip of the Appendix.) Clearly, the bisepara- 
bility condition of p7al) is the same as the one of p2ap of 
the criterion on spin-observables but condition of p7bp 
is strictly stronger than the one of (j32bp . (On the 5 = 
noisy W state it gives bound 9/17.) The full-separability 
condition of (|38ap is the same as the one of (|33ap of the 
criterion on spin-observables (and the one of (I17ap of par- 
tial transposition criterion as well) but the condition of 
(I38cp is weaker than the one of (I33bp of the criterion on 
spin-observables. Hence at this point these criteria are 
stronger for biseparability but weaker for full separability 
than the criteria on spin-observables for our state. But 
we have another full-separability condition: (j38bp can be 
stronger in a region than the ones based on the partial 
transposition criterion. The states of parameters in this 
region are entangled ones of positive partial transpose, 
no pure state entanglement can be distilled from them. 
The borders of the domains in which these conditions 
hold and the region of PPTESs can be seen in Fig. HI 

Now we show a representing matrix of the region of 
PPTESs determined by p7a|) . (fTZb)) and (l38b)) . One can 
check that the state of parameters g = 1/5, w = 1/5 is 
contained by this set and the explicit form of Q for this 
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point is 



can detect 2?™^ and 



PPTES 



120 



21 



17 8 • 8 

8 17 • 8 

• • 9 • 

8 8 • 17 



12 



12 



21 



(39) 



VI. SLOCC CLASSES OF GENUINE 
THREE-QUBIT ENTANGLEMENT 

A fully entangled three-qubit pure state can be ei- 
ther of GHZ-type or of W-type in the sense of 
Stochastic Local Operations and Classical Communica- 
tions (SLOCC, [l3|): vectors of these two different types 
can not be transformed into each other by local invert- 
ible operations. These fully entangled vectors can 
be classified by the so called three-tangle ti23(V') UM 
'''123 ("0) 7^ exactly for the GHZ-type vectors. 

In Ref. [15] Acin et. al. have investigated the classifi- 
cation of mixed three-qubit states in the sense of SLOCC 
and they have shown that Class 1 of fully entangled states 
can be divided into two subsets, namely the ones of GHZ 
and W-type entanglement, by the following definitions. 
A state is of W-type (2?^) if it can be expressed as a 
mixture of projectors onto 2-separable and W-type vec- 
tors (therefore is also a convex set) and GHZ-type 
vector is required for a GHZ-type mixed state. Hence 
the following holds: 



GHZ 
3 



pi -sop _ p^ 



(40) 



Let Class W the set Vf 



\ X>3""°P and Class GHZ the 
set Vf^^\Vf, so Class 1 = Class WU Class GHZ. The 
mixed state extension of T123 



Ti23(£') = mill 



PjTi23(V'i) 



E 



Pi\ipi){tpi 



(41) 

is a good indicator for Class GHZ: Ti23{g) ^ exactly 
for Class GHZ. 

A method to determine to which class a given mixed 
state belongs is the use of witness operators. Q A her- 
mitian operator W is a witness operator for a convex 
compact set C C I? AT if < TrW^cr for all u e C and 
there exists g ^ C for which Ti Wg < 0. Hence the neg- 
ativity of expectation value of the observable W bears 



witness that the state does not belong to C. In Ref. [12 1 



there have been given some witnesses for Vj^ and V^^^ : 



W, 



GHZ 



I GHZ) (GHZ I 



(42) 



Ww^ = f-m{W\, 



W2 



|GHZ)(GHZ| 



(43a) 
(43b) 



can detect V^. With these we have 



gev^ 



< TriyGHz£'= (20rf- 2g + 9w)/4 
= (5 - 7g + w)/8 



(44) 



and if the inequality is violated then g G Class GHZ, as 
well as 



^ ^2— Sep 



< TrVKwi(? = il3d + Ag-3w)/3 
= (13-t-3g-21w)/24, 

< TtWw^Q = {6d-2g + 3w)/2 
= {3-7g + w)/8 



(45a) 
(45b) 



and if either or both of the inequalities is violated then 
g e Class 1. In Fig. IHlwe plot the lines on which these 
inequalities are saturated. It can be checked that (|45ap 
and (|45bp gives weaker condition for biseparability than 
P7ap and (|37b[) of the previous Section. We can conclude 
that all the states in the blue/grey domain belong to 
Class GHZ, and the biseparable states are enclosed by the 
blue/black curves, however, both type of fully entangled 
states can be here too. 

The equality in (j44p gives an "upper bound" for the 
border of Class GHZ. (See blue/grey domain in Fig. [9]) 
Fortunately, we have a possibility to give also a "lower 
bound" for that, thanks to the results of Lohmayer 
et. al. [19]. They have studied the GHZ-W mixture {d = 
0) and they have found that there exists a decomposition 
of projectors onto vectors of vanishing three-tangle if and 
only if < 5 < ffo = 4 • 2^^^/{S + 4 • 2^^^) = 0.626851 . . . , 
hence for these parameters the mixed state extension ([^T|) 
of the three-tangle is zero. If we mix the states of this 
interval with white noise then the three-tangle remains 
zero and neither of these states can belong to Class GHZ. 
So we can state that 



g e 



w < 



4-21/3 



9- 



(46) 



which holds under the green/light grey line of Fig.|9l This 
condition is quite weak, but we can make it stronger. 
Recall that on the w = line (noisy GHZ state) g e 
Class 1 if and only if 3/7 < g < 1. (See Section IHIl) 
So the convexity of restricts Class GHZ to be inside 
the triangle defined by the vertices {g = 3/7, w = 0), 
(g — l,w — 0) and {g = gQ,w ~ 1 — go)- (Union of tinted 
domains in Fig.|9l) So we can conclude that all the states 
in the blue/grey domain belong to Class GHZ, and the 
border of Class GHZ is in the red/light grey domain of 
Fig.\E 
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FIG. 9. (Color online) Criteria on three-partite entanglement 
classes for state Q on the p-w-plane. (Blue/black straight 
lines: the borders of domains inside Eqs. (|45a|l - (|45b|) hold, 
blue/black curves of second order: the borders of domains 
inside Eqs. (I37a)l - (f37bll hold, copied from Fig. \8\ Eq. ^ 
holds inside the blue/grey domain, and the border of Class 
GHZ is inside the red/light grey domain. The inequalities 
hold on the side of the curves containing the origin. Eq. (|46p 
holds under the green/light grey line.) 



VII. CONCLUSIONS 

In this paper we have investigated the noisy GHZ-W 
mixture and demonstrated some necessary but not suf- 
ficient criteria for different classes of separability. With 
these criteria we can restrict these classes into some do- 
mains of the 2-dimension simplex. It has turned out that 
the strongest conditions was (|17a|) . (j38b| and (|17bp for 
full separability, (|17ap and (|17b|) for the union of Classes 
2.8 and 3 and (I37ap and (j37bp for biseparability. These 
have been obtained from the partial transposition crite- 
rion of Peres [l| and the criteria of Giihne and Seevinck 
01 dealing with matrix elements. Only these latter crite- 
ria have turned out to be strong enough to reveal a set of 
entangled states of positive partial transpose. (The set 
of these states can be given by the conditions of (I17ap . 
(|17bD and (|38bp . An example is given in Eq. (15^1).) 

Besides this, some remarkable coincidences have also 
appeared: some parts of some bipartite separability cri- 
teria have proved to be necessary and sufficient for sepa- 
rability classes of the GHZ-white noise mixture, (e.g. the 
majorisation criterion and the entropy criterion in the 
a — oo limit as well, and reduction criterion.) This 
is interesting because e.g. the majorisation criterion — as 
our knowledge — does not state anything about a density 
matrix which is majorized by only either of its subsys- 
tems. We do not think that this would be more than a 
coincidence, however, the GHZ state is very special so it 
is an interesting question whether this can be generalized 
to the n-qubit noisy GHZ state. 



Another interesting observation was that the two set- 
tings/measurement vectors strong in detection of GHZ 
and W state are related by local unitary Hadamard trans- 
formation in the criteria on spin-observables of Seevinck 
and Uffink (Section IV Bp and also in the criteria of 
Gabriel et. al. (Section lV CI) . The transformation on set- 
tings/measurement vectors can also be written on the 
state Q I—)- {H'^^y qH®^ , which means that we can use 
the same measurements on the transformed state for the 
detection of W state as for the detection of GHZ state. 

We have also investigated the SLOCC classes of fully 
entangled states and we have given restrictions for Class 
GHZ. 
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Appendix A: Matrices 

In this Appendix we show the density matrix g given 
in Eq. (jj]), its partial transpose, its reshufflings, and its 
marginals explicitly. Note that we use the renormalized 
parameters d = d/S, g = g/2, w = w/3 and the con- 
straint d + g + w = 1 holds. 



d + g 



d + w w ■ w 
w d + w ■ w 
d ■ 

w w ■ d + w ■ 
d 

d 



d + g_ 



, (Al) 



d+g ■ ■ ■ ■ WW 

d + w w 
w d + w ■ 

d g 
g d + w ■ 

w ■ ■ ■ ■ d 

w ■ ■ ■ ■ ■ d 



d + l 



„23 



2d + g + w 



2d + w w 
w 2d + w 



2d + g 



(A2) 



(A3) 
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4:d + g + 2w 



Ad + g + w 



(A4) 



d + w 



d ■ ■ —w 
d ■ —w 

■ d + g _ • 

-w —w ■ d + g 



d + w w 
w d + w 



-9 



(A5) 



3d+2w 



3d + g + w —w 
—w 3d + g + w 



-w 
-w 



3d + g + 2w 



-w 



-w 



3d + g 



3d + g + w 



3d + g + w 



3d + w 



(A6) 



R{q) 



d + g ■ ■ ■ ■ d + w w ■ ■ w d + w ■ ■ ■ ■ d 

■ ■ g w ■ ■ ■ w 

W W ■ ■ • .... . . g . . 

d + w ■ ■ ■ ■ d .... d ■ ■ ■ ■ d + g 



(A7) 



R'{g) = 



d + g ■ - d + w 

■ w ■ ■ g 

w ■ ■ w 

d + w ■ ■ d 

w ■ ■ d+w 

w ■ ■ ■ 



w 



d + g 



(A8) 
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Appendix B: Permutation-invariant states of Cleisses invariant 2-separable one: 
2.1 and 2.8 



To see that Class 2.1 is not empty for permutation- 
invariant states in general we will show an explicit ex- 
ample. Let |/9o) = (|00) + \11))/V2 be the Bell-state, 
then the uniform mixture of the rank one projectors to 

the subspaccs |0)i (g) |/3o)23, |0)2® |/3o)3i and |0)3 (g |/3o)i2 
gives a state which is by construction a permutation- 



^Class 2.1 



1-11 



(Bl) 
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It can be easily checked that its partial transpose is not 
positive, so it is not a2-separable, (the partial transpo- 
sition criterion see in Section IIV C[) hence it is in Class 
2.1. 

An example for a permutation-invariant state in Class 
2.8 is given in Eq. (14) of with a = b = ^: 



Clearly, the last eigenvalue is the smallest one. If 4((i -|- 
w)2 < {2d+g){2d+g+w) then Af ^ hence C{g'^^) = 0. 
If 4(J -t- w)2 > (2d + g){2d + g + w^, 
C{g^^) can be nonzero. It turns out that 



then Aj 



Ag and 



9 



Class 2.8 



2 + 3(a+i) 



1 



1 



1 



1 



(B2) 

where < a. This state is entangled if and only if a ^ 1, 
and 2.8 g 25^2 g| 



C(£»23) = 2w - 2y^(2J+.g)(2j+.g-hw) 



3" 



2^/3 



v/(H-ff-u.)(3 + 35 + u;) 



(C3) 



if < - {2d + g){2d + g + w) = (-9g2 -t- lOw^ + 6gui - 
18g -f 6w - 9)/122, otherwise C{g^^) = 0. (Fig. [TO]) It 
takes its maximum 2/3 in g = 0, w = 1, i.e. in pure W- 
state. For the GHZ-W mixture {d = 0) we get back the 
result of Ref. [l9,]. 



Appendix C: Wootters concurrence 



For mixed states the only measure of entanglement 
known explicitly is the Wootters- concurrence of two 
qubits: [11 [11] 

C(w) = max{0, A]; - A^ - A^ - A:^}, (CI) 

where w is a two-qubit density matrix and A^s are the 
eigenvalues of the positive matrix -y/ ^Ju)Cj^/u) in non- 
increasing order. Here uj denotes the spin- flipped density 
matrix: Cj = [a2 ® <72)uj*{(J2 ® cr2) where f72 = [? V] 
and UJ* denotes the complex conjugation. Equivalently 
AiS can be calculated as the square roots of the positive 
eigenvalues of the matrix cDw. 

Let us calculate now the concurrence given in Eq. (jCip 
of the matrix g^^. (See in Eq. (jA3[) of the Appendix.) 
Since the spin-flip in the two-qubit case means transpose 
with respect to the antidiagonal then multiplication of 
neither diagonal nor antidiagonal entries by —1, one can 
easily get: 

Spect(^23(?23) ^ 

{ 4(d 4-^)2 =4(3- 3.g + 5w)V24^ (C2) 
{2d->rg){2d + g + w) 

= 12(1 +g- w)(3 + 3g + 'u;)/24^ 
{2d + g){2d + g + w) 

= 12(1 + .9 - w){3 + 3g + w)/2A^, 
4(P =36(1-5-^)2/24^ } 




[1] A. Peres, |Phys. Rev. Lett.| 77, 1413 (Aug 1996). 



[2] O. Giihne and M. Seevinck, New Journal of Physics 12, 



18 



053002 (2010). 

R. Horodecki, P. Horodecki, M. Horodecki, and [20 
K. Horodecki, [RevTModrP^s] 81, 865 (Jun 2009). 
M. A. Nielsen and I. L. Chuang, Quantum Computation [21 
and Quantum Information, 1st ed. (Cambridge Univer- 
sity Press, 2000) ISBN 0521635039. [22 
C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, 
A. Peres, and W. K. Wootters, |Phys. Rev. Lett^l TO, 1895 [23 
(Mar 1993). 

0. Giihne and G. Toth, [Physics Reports] 474, 1 (2009). [24 
R. F. Werner, Phys. Rev. A" 40, 4277 (Oct 1989). 
M. Horodecki, P. Horodecki, and R. Horodecki, [25 
Physics Letters A 223, 1 (1996), ISSN 0375-9601. 

1. Bengtsson and K. Zyczkowski, Geometry of Quan- [26 
turn States: An Introduction to Quantum Entanglement 
(Cambridge University Press, New York, NY, USA, 2006) [27; 
ISBN 0521814510. 

M. Seevinck and J. Uffink, |Phys. Rev. A| 78, 032101 (Sep [28 
2008). 

C. H. Bennett, D. P. DiVince nzo, T. Mor, P. W . Shor, [29 
J. A. Smolin, and B. M. Terhal, [Ply8. Rev. LettT] 82, 5385 
(Jun 1999). 

A. Acm, D. Brufi, M. Lewenstein, and A. Sanpera, [30 
Phys. Rev. Lett. 87, 040401 (Jul 2001). 

C. H. Bennett, S. Popescu, D. Rohrlich, J. A. Smolin, and [31 
A. V. Thapliyal, Phys. Rev. A 63, 01 2307 (Dec 2000 ). 
W. Diir, G. Vidal, and J. I. Cirac, |Phys. Rev. A] 62, [32 
062314 (Nov 2000). 

V. Coffman, J. Kundu, and W. K. Wootters, 
Phys. Rev. A 61, 052306 (Apr 2000). [33 
A. Uhlmann, Phys. Rev. A 62, 032307 (Aug 20 00). 
S. Hill and W. K. Wootters, [Phj^Rev. Lett.] 78, 5022 [34 
(Jun 1997). 

W. K. Wootters, [Phys. Rev. LetT !l80, 2245 (Mar 1998). [35 
R. Lohmayer, A. Osterloh, J. Siewert, and A. Uhlmann, 



'Phys. Rev. Lett.' 97, 260502 (Dec 2006). 

P. Horodecki, Physics Letters A 232, 333 (1997), ISSN 

0375-9601. 

W. Diir, J. I. Cirac, and R. Tarrach, [Fhys. Rev. LetE] 83, 
3562 (Oct 1999). 

M. A. Nielsen and J. Kempe, |Phys. Rev. Lett. [ 86, 5184 

(May 2001). 

R. Horodecki and M. Horodecki, [Phys. Rev."A| 54, 1838 
(Sep 1996). 

R. Horodecki, P. Horodecki, and M. Horodecki, 
Physics Letters A 210, 377 (1996), ISSN 0375-9601. 

B. M. Terhal, Theoretical Computer Science 287, 313 
(2002), ISSN 0304-3975. 

K. G. H. Vollbrecht an d M. M. Wolf, 
[Journal of Mathematical Physics[ 43, 4299 (2002). 
M. Horodecki, P. Horodecki, and R. Horodecki, 
[Ph y s. Rev. Lett Ti 80, 5239 (Jun 1998). 

C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. 
Wootters, Phys. Rev. A 54, 3824 (Nov 1996). 

L. Clarisse, "Entanglement distillation; a discourse on 
bound entanglement in quantum information theory," 
(Dec 2006), | ^Xi?iqu ant-ph/061207 2?T| 
M. Horodecki and P. Horodecki, P hys. Rev. A, 59, 4206 

(Jun 1999). 

N. J. Cerf, C. Adami, and R. M. Gingrich, [Fhys. Rev. A[ 
60, 898 (Aug 1999). 

M. Horodecki, P. Horodecki, and R. Horodecki, Open 
Systems & Information Dynamics 13, 103 (2006), ISSN 
1230-1612, 10.1007/sll080-006- 7271-8. 
L. Clarisse and P. Vocjan, Quantum Information and 
Computation 6, 277 (May 2005). 

A. Gabriel, B. C. Hiesmayr, and M. Huber, Quantum 
Information and Computation 10, 829 (2010). 
M. Huber, F. Mintert, A. Gabriel, and B. C. Hiesmayr, 
|Phys. Rev. Lett. 104, 210501 (May 2010). 



